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Abstract. Recently, A. I. Aptekarev and his collaborators found a sequence of ratio- 
nal approximations to Euler's constant 7 denned by a third-order homogeneous linear 
recurrence. In this paper, we give a new interpretation of Aptekarev's approximations 
to 7 in terms of Meijer G-functions and hypergeometric-type series. This approach 
allows us to describe a very general construction giving linear forms in 1 and 7 with 
rational coefficients. Using this construction we find new rational approximations to 7 
■ generated by a second-order inhomogeneous linear recurrence with polynomial coeffi- 

' cients. This leads to a continued fraction (though not a simple continued fraction) for 

Euler's constant. It seems to be the first non-trivial continued fraction expansion con- 
' vergent to Euler's constant sub-exponentially, the elements of which can be expressed 

as a general pattern. It is interesting to note that the same homogeneous recurrence 
generates a continued fraction for the Euler-Gompertz constant found by Stieltjes in 
1895. 

> 

o . 

(N . 1. Introduction 

In 1978, R. Apery [3, 31] stunned the mathematical world with a proof that £(3) = 
Y^Li ^3 is irrational. Since then many different proofs of this fact have appeared in 
the literature (see [15] and the references given there). The main idea of all the known 
proofs is essentially the same and consists in constructing a sequence of linear forms 
I n = u n ((3) — v n , n = 0, 1, 2, . . . , which satisfy the following conditions: 

>< : limsup \I n \ 1/n < (V2 - l) 4 = 0,0294372. 

5-H 
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/„ ^ for infinitely many n, and u n G Z, 2Df l v n G Z, where D n is the least common 
multiple of the numbers 1,2, ... ,n. If we suppose that ((3) is a rational number a/b, 
then 2bDf l I n is a non-zero integer for infinitely many n and, on the other hand, it tends 
to zero as n — >■ 00 (since Dl/ n — > e and e 3 (y/2 — l) 4 = 0.591 ... < 1), which is a 
contradiction. 

The diversity of all the proposed proofs of the irrationality of ((3) is presented by 
different interpretations and ways of obtaining the linear forms I n and its various rep- 
resentations. Apery showed that the sequences u n and v n are given by the formulas 
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and satisfy the second-order linear recurrence relation 

(n + l)V+i - (34n 3 + 51n 2 + 27n + 5)y n + n 3 y n ^ = 

with the initial conditions u = l,Ui — 5, v — 0, v\ = 6. This implies immediately that 
^ is the n-th convergent of the following continued fraction: 

_ 6 1 64 n 6 

^ ^ ~ 5- 117-535 - - 34n 3 + 51n 2 + 27n + 5 - "" 

The shorter proof of Apery's theorem has been found in 1979 by F. Beukers [8], who 
used multiple Euler-type integrals and Legendre polynomials. Beukers discovered the 
following elegant representation for the linear forms I n : 

(1) 

P n (x)P n (y) \og(xy) fff x n (l - x) n y n (l - y) n w n {\ - w) n 



2I " = J J — —i — ixdv = JJJ (i-(i -*,)»)-« dxdydw ' 

[0,1] 2 [0,1] 3 

where P n (z) = (z n (l — z) n )^ is the Legendre polynomial. 

In 1996 inspired by the works of F. Beukers [9] and L. Gutnik [16], Yu. V. Nesterenko 
[28] proposed another proof of the irrationality of £(3) and a new expansion of this 
number into continued fraction. His proof was based on the hypergeometric type series 

P) " 5 £ = "5 g S P(«- B gg +B + i) L = U ' m ~ 

that can be written (by the residue theorem) as a complex integral or a Meijer G- 
function (see [27, §5.2], for definition) 

c+ioo c+ioo 

1 f n „M \ 2 , i f T 2 (n + 1 - s)r 4 (s) , 

as 



(3) fc=1 



T 2 (n+ 1 +s) 



C—ICO 



_ r 4,2 / n + 1, n + 1 

~ 4,4 V 0, 0, 0, 



here c is an arbitrary real number satisfying < c < n + 1 and -R n (t) is a rational 
function defined by 

7? ^ (*-l) 2 ---(*-n) 2 

Using an identity of Bailey (see [7, Section 6.3, formula (2)]) it can be shown (see [45]) 
that the complex integral in (3) is equal to the very-well-poised hypergeometric series 
(see [15, §2], for definition) 

fA . n! 7 (3n + 2)! / 3n + 2, \n + 2, n + 1, . . . , n + 1 

W 1^ I 1 \ 15 7^6 3 



(2n+l)! 5 U V |n+ l,2n + 2,...,2n + 2 

This series was first introduced by K. Ball in attempts to give an elementary proof of 
Apery's theorem similar to the proof of the irrationality of the number e by Fourier. As 
we can see now, the hypergeometric construction (2)-(4) appeared to be more trans- 
parent for generalizations. A generalization of the series (4) allowed T. Rivoal [33] to 
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prove his remarkable result on infiniteness of irrational numbers in the set of odd zeta 
values C(2n + 1), n > 1. The best finite version of it belongs to W. Zudilin [44] and 
states that at least one of the four numbers £(5), £(7), £(9), £(11) must be irrational. 
Another generalization of (2), (3) to values of the polylogarithmic function 



enabled the second author [24, 25] to prove some results on arithmetical nature of 
values of the polylogarithm at rational points lying inside the disk of convergence far 
enough from zero. Namely, for any integer q ^ 0, 2 at least one of the two numbers 
L 2 (l/q), L 3 (l/q) is irrational and for any integer q G (— oo, —4] U [7, +oo) at least one 
of the three numbers L 3 (l/g), L 4 (l/g), L 5 (l/q) is irrational. For small \l/q\ we have 
more precise results due to M. Hata [18, 17] and G. Rhin and C. Viola [32]: L 2 (l/g) is 
irrational for any q G (— oo, —5] U [6, +oo), and L 3 (l/q) is irrational for any integer q 
with \q\ > 1038. 

Further extension of the hypergeometric approach led to the discoveries of second- 
order linear recursions and Apery-like continued fractions for other remarkable number- 
theoretic constants including Catalan's constant [46], 7r 4 , [11, 47, 39], 7rcotli7r [20]. For 
other related examples we refer the reader to [2, 43, 48]. 

Euler's constant was first introduced by Leonhard Euler in 1734 as 



It can be considered as an analogue of the value "C(l)" of Riemann's zeta function if 
we compensate the partial sums of the divergent harmonic series by the natural log- 
arithm. It is not known whether 7 is an irrational or transcendental number. The 
question of its irrationality remains a famous unresolved problem in the theory of num- 
bers. Even obtaining good rational approximations to it was unknown until recently. 
First such approximations defined by a third-order linear recurrence were proposed by 
A. I. Aptekarev and his collaborators [4] in 2007. More precisely, the numerators p n and 
denominators q n of the approximations are positive integers generated by the following 
recurrence relation: 

(16n - 15)y n+ i = (128n 3 + 40n 2 - 82n - Ah)y n 

- n 2 (256n 3 - 240n 2 + 64n - 7)y n - 1 + n 2 (n - l) 2 (16n + l)y n _ 2 

with the initial conditions 




oo 





- log n) = 0.57721566490153286 .... 



Po = 0, 
<?o = 1, 

and having the following asymptotics: 



Pi = 2, 

qi = 3, 



P2 = 31, 

h = 50 
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The remainder of the above approximations is given by the integral [5] 

PCX) 

(7) / Q n (x)e~ x log(x) dx = p n -79™, 

Jo 

where 

Qn(x) = \-^—(x n {x n (l - x )2n+l e -)W)W 

nr 1 — x 

is a multiple Jacobi-Laguerre orthogonal polynomial on [0, 1] and [1, +00) with respect 
to the two weight functions W\(x) = (1 — x)e~ x , W2{x) = (1 — x) log(x)e~ x . The integral 

(7) can also be written as a multiple integral (see [22, Lemma 4]) 

poo poo poo n n( _ i \2n+l p -x 

(8) / Q n (x)e- x log(x)dx= / / y -hi e dxdy , 

KJ Jo Jo Jo (xy + i) n+1 (y + i) n+1 

The integrality of the sequences p n and q n is not evident and can not be deduced directly 
from the recurrence equation (5). Tulyakov [41] proved independently that p n and q n 
are integers, by considering a more "dense" sequence of rational approximations to 7. 
The authors [23] found explicit representations for p n and q n : 

n , x 2 n , x. 2 

(9) Qn = J2[ n k ) ( n + k V-i Pn = J2[ n k ) (n + k)\(H n+ k + 2H n _ k -2H k ), 

k=0 ^ ' k=0 ^ ' 

here H n = X)fc=i | * s the n_ th harmonic number, if := 0. Formulas (9) imply that q n 
and p n are integers divisible by n! and respectively. Although the coefficients of the 
linear forms (6) can be canceled out by the big common factor it is still not enough 
to prove the irrationality of 7, since the linear 7-forms with integer coefficients: 

PnD n q n D n 

: r 7 G Z + Z7 

do not tend to zero as n tends to infinity: 



PnDn _ q n Dn _ Q^n n n-l/4 e -V2n~^ 



nl n\ 

Nevertheless, the sequence j 1 provides good rational approximations to Euler's constant 
^-7 = 27re- 2 ^(l + 0(n- 1 / 2 )) as n ->• 00. 

Qn 

In 2009, T. Rivoal [35] found another example of rational approximations to the Euler 
constant 7, by using multiple Laguerre polynomials 

1 

n! 2 

His construction is based on the following third-order recurrence: 



M*) = z ^e x (x n (x n e- x Y n) Y n) . 



{n + 3) 2 (8n + ll)(8n + 19)y n+3 = (n + 3)(8ra + ll)(24n 2 + 145n + 215)y n+2 
- (8ra + 27)(24n 3 + 105n 2 + 124n + 25)y n+1 + (n + 2) 2 (8n + 19) (8ra + 27)y n , 



which provides two sequences of rational numbers P n and Q n , n > 0, with the initial 
values 

P = -1, P l = 4, P 2 = 77/4, 

Qo = 1, Q x = 7, Q 2 = 65/2 

such that 7f- converges to 7. The sequences P n , Q n satisfy the inclusions (see [22, 
Corollary 5]) 

n\ Q n , n\ D n P n e Z 
and provide better approximations to 7 

< Coe -9/2nV3 + 3/2nV3 ^ | Qf| | = O^'""- 1 ^ &g „ ^ 

Unfortunately, this convergence is not fast enough to imply the irrationality of 7. 

In this paper, we give a new interpretation of Aptekarev's approximations to 7 in 
terms of Meijer G-functions and hypergeometric-type series, which can be considered 
as an analog of the complex integral (3) and series (2) for Euler's constant. This 
approach allows us to describe a very general construction giving linear forms in 1 and 7 
with rational coefficients. Using this construction we find new rational approximations 
to 7 generated by a second-order inhomogeneous linear recurrence with polynomial 
coefficients. This leads to a continued fraction (though not a simple continued fraction) 
for Euler's constant. It seems to be the first non-trivial continued fraction expansion 
convergent to Euler's constant sub-exponentially, the elements of which can be expressed 
as a general pattern. It is interesting to note that the same homogeneous second-order 
linear recurrence generates a continued fraction for the Euler-Gompertz constant found 
by Stieltjes in 1895. 



p 

1 n 
Qn 



2. Analogs of hypergeometric type series and complex integrals for 

Euler's constant 

Note that the first attempt to generalize series (2) to find suitable approximations 
for Euler's constant 7 was made by Sondow [38]. He introduced the following series: 

and proved that 

(10) I n = 

where 

n k—1 / \ 2 

Ln = 2j2J2{ n t ) (Hn-i-H^login + k). 

k=l i=0 ^ ' 

Prevost [30] showed that the construction (10) is a consequence of the Pade approxi- 
mation to the function (log it) /(l — u) and proposed simpler approximations to Euler's 
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- E 



u=n+l 



v \x(x + 1) • • • (x + n) 



dx 



4 "n _1 ! 



constant given by 



J n 



U 



l Pn{u) 



1 1 
+ " 



logw 



u 



du = ^ - 2H n + L n>m , 



where 



is a linear form in logarithms and P n {u) is the Legendre polynomial. Rivoal [34] men- 
tioned on the following connection 



7 



S n — 



u 



1 1 

+ 



log u 1 — u 



where S n 
shown in [19] that 



Y2=i \ - lo g( n + 1)) and obviously, 7 = lim T1 



S n . In particular, it was 



Jn,m 7 



E(-D" + '(I) 



n + k 
k 



~>k+n—m- 



This gave a new view on the approximating formulas for Euler' constant 7 found much 
earlier by different methods by C. Eisner [12]: 



7 



k=o v 7 



fc + n + r — 1 
k + T-1 



k+T-l 



< 



2-(T)' 



TEN, 



and led to proving many similar formulas for 7 and some other constants (see [13, 
19, 21, 34]). Unfortunately, this type of approximations involving S n contains linear 
forms in logarithms and this enables only to obtain conditional irrationality criteria for 
7, accelerate convergence of the sequence S n and provide ways for fast calculation of 
Euler's constant. 

In this section, we obtain new representations for Aptekarev's linear form p n — 
distinct from (7) and (8) which can be considered as analogs of the hypergeometric-type 
series (2) and complex integral (3). It can be easily done by using explicit formulas (9). 

Proposition 1. For each n = 0, 1, 2, . . . , the following equality holds: 



fn ■ Pr, 



d 



k=0 



r(n + t + i) 



dt \T 2 (t + l)T 2 (n-t + l) 



t=k 



Proof. The straightforward verification shows that 



d 



'ID 



T(n + t + r 



T(n + t 



dt \T 2 (t + l)r 2 (n -t+l)J T 2 (t + l)T 2 (n - t + 1) 

x + * + !)- 2ip(t + 1) + 2^(n - * + 1)) , 



where ip(z) 



EM. 

r(z) 



is the logarithmic derivative of the gamma function, also known 



as the digamma function. Summing (11) for t = 0,1, ... ,n and using the well-known 
properties of the digamma function 

^(l) = -7, ^(n + 1) = H n - 7, n>l, 
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we have 



E 

fc=0 



d 



T(n + t + l) 



dt \T 2 (t+ l)r 2 (n-t + 1) 



t=k 



E 

k=0 



T(n + k+ 1) 



P(A; + l)r 2 (n-A; + l) 



(ip(n + k + i) 



: + 1) + 2^(n - k + 1)) 



n / \ 2 

n 



E 



fc 



n + k)\{H n+k - 2H k + 2H n _ k - 7) = p n - 7^, 



fc=0 

as required. 

Proposition 2. For eac/i n = 0, 1, 2, . . . , we have 



□ 



(12) 



c+ico 



I> + t + l) 



T 2 (t + l)T 2 (n-t + l) 
n] 2 rc+i°°r(n + t + l)r 2 (t-n) 



1) Vsin7rt/ 



nr 



q n = 



r 2 (t + i) 

Yin + t + 1) e i7Tt 



r, ,2^0,3 f-n,n + l,n+ 1 



r 3,2 



0, 



7T 



2vri 7 L T 2 (t + l)r 2 (n - t + 1) sin vrt 



dt 



(13) 



^-l)"n! 2 f T(n + t+ l)r(t-n) 



„l7rt 



2?ri y L r 2 (t + i)r(n-t + 1) 

: -l)»nPG^ (-"'"+ 1 , " +1 |-l), 



dt 



where c is an arbitrary real number satisfying c > n and L is a loop beginning and 
ending at —00 and encircling the points n, n — 1, n — 2, . . . exactly once in the positive 
direction. 

Proof. We first note that the equality of two complex integrals in (12) and (13) follows 
easily by the reflection formula for the gamma function: 



r(t - n)T(n - t + 1) 



(-iyv 

shi7rt 



The third equality in both formulas (12) and (13) follows by the definition of the Meijer 
G-functiom To prove the first equality in (12), we consider the integrands of the complex 
integrals (12) on the rectangle contour with vertices c ± iN, —N — 1/2 ± iN, where N 
is a sufficiently large integer, N > c. Then, by the residue theorem, we have that the 
integral 



(14) 



2ni 



c+iN 



/ + / + / 

Jc+iN J-N-\+iN J-N-\ 



AT-i-iJV 



rc-iN 



+ 



c-iN 



N-±+iN 
7 



-iN , 



T(n + t+ l)r 2 (t-n) 

r 2 (t + i) 



dt 



is equal to the sum of the residues of the integrand at integer points t = k, < k < n. 
Using the expansion 

7T \2 1 

+ 0(1) 



^sin7rt/ (t — k) 2 

in a neighborhood of the integer point t = k we obtain that the integral (14) is equal to 
(15) 

T(n + t+ 1) 



V res ( T ^ + t + ^ (^-Y) =T±( 

^t=k 1 T 2 (t+ l)r 2 (n-t + 1) Vsinvrt; J ^ dt V 

fc=0 \ v 1 / V J / fc=0 \ 



r 2 (t + i)r 2 (n-t + i) 



t=k 



Since on the sides [c + iJV, -iV - 1/2 + i/V], [-iV - 1/2 + iN, -N - 1/2 - iN], [-N - 
1/2 — iN, c — iN] of the rectangle we have \t\ — 0(N), it follows that 



(16) 



T 2 {t-n) 



r 2 (t + i) 



o 



\t 2 (t- I) 2 ■ ■ ■ (t - n) 2 \ \N 2n + 2 



For large \z\ the asymptotic expansion of the gamma function is [27, §2.11] 



(17) 



logr(z) =( z -^\logz-z + ± log(27r) + 0(\z\- 1 ), 



where | argz| < ti — e, e > and the constant in O is independent of z. Then for 

t — x ± iiV, — iV < rr < c, we have 

(18) 

|r(n+t+l)| = \T(x+n+l±tN)\ = 0(e^ +n+ ^ losN ^ Nars{x+n+1±iN) ) < 0{N c+n+ ^e~^ N ). 

On the segment \—N — 1/2 + iN, — iV — 1/2 — iN] we use the trivial estimate 
(19) 

|r(n+f+l)| < |r(Re(n+t+l))| = |r(n+l/2-iV)| = v{N + \ /2 _ n) = 0(e' N ^ N+N ). 

Summarizing (16), (18), (19) and letting iV tend to infinity in (14), by (15) and Propo- 
sition 1, we get the first equality in (12). 

Now prove the first equality in (13). For this purpose, suppose that C\ and C 2 are 
points of intersections of the loop L with the vertical line Re (t) = —N — 1/2 and 
consider a closed contour L* oriented in the positive direction and consisting of the 
segment CiC 2 and the right part of the loop L connecting the points C\ and C 2 , which 
we denote by CiC 2 - Then, by the residue theorem, we obtain 
(20) 

n\ 2 [ T(n + t + 1) e i7Tt dt _ |r A / T(n + t + 1) e i7Tt vr \__ 
~2i J L , Y 2 (t + l)r 2 (n - t + 1) imTt ~ f^t=t \T 2 (t + 1)T 2 (n - t + 1) irn^rt ) ~ 

On the other hand, we have 

(21) / = /_+/ • 

J L* J C\Ci JC2C1 

Since 

F(n + t + 1) e int ti 1 T 2 (t-n) , ,. 9 . f 

5: ! ! L — 5: -Yin 4- / 4- \\{e — 1) 

r 2 (t + l)r 2 (n-t+l) shiTrf 2m F 2 (t + 1) V n ; 



and the function e 2mt — 1 is bounded on the vertical segment CiC 2 , by (16), (19) and 
(21) we get 



[ = [ +0{e - mog N + Ny 

Jl* Jcu% 



Now letting N tend to infinity and taking into account (20), we get the required formula. 

□ 

Remark. Note that the asymptotics of linear forms p n — ^q-n and denominators q n can 
be also obtained by application of the saddle-point method to the complex integrals 
(12) and (13). 

More generally, let us consider an arbitrary function F(n, t) of the form 

(22) F(n,t) = tX 3 r 7> 

V 1 V ' n;=i r(c,-n + djt + 1) 

where cij, bj, Cj, dj G Z, J^=i ^ Yl]=i dj, anc ^ a ^ gamma values in the numerator are 
well defined for < t < M{n). We say that Y{an + bt + l) is well defined if an + bt + 1 
is not a negative integer or zero. 



Proposition 3. Let F(n,t) be defined as above, M(n) be a non-negative integer, and 

M(n) 
t=0 

Then for each n — 0, 1, 2, . . . , we have F n = p n — ^q n with 

s u M{n) 

3=1 j=l k=0 

and 

M(n) s u 

k=0 3=1 3=1 

Proof. Differentiating F(n,t) with respect to t and summing over t — 0, 1, . . . ,M(n), 
we have 



t=0 k=0 11 J 

X 



bjip{ajn + bjk + 1) — djip{cjn + djA; + 1) j 



nLi(«^ + bjk)\ / * " x 

= 2^ n" (V^TTfrV 12-^ b 3 {H ajn+bjk - 7) - 2^ dj(H Cjn+djk - 7) J = p n - 7 g n . 

□ 



3. A SECOND-ORDER INHOMOGENEOUS LINEAR RECURRENCE FOR EULER'S 

CONSTANT 

In this section, we consider application of Proposition 3 to the function 

F(lM) = r( t + S(t- t + i) - »ez,»>a 



Then for 



dt 

t=o 

we have F n = p n — jq n with 

n / \ 2 n / \ 2 

(23) 9n = Xl(j A;! ' P» = Z)( J k\(2H n _ k -H k ), n = 0,l,2,.... 

fc=0 ^ ' fc=0 ^ ' 

Lemma 1. T/ie sequence {g n }^l is a solution of the second-order homogeneous linear 
recurrence 

(24) q n+2 - 2(n + 2)q n+l + (n + l) 2 g n = 

with the initial values q — 1, qi — 2, and £/ie sequence {p n }^L is a solution of the 
second-order inhomogeneous linear recurrence 

n 

(25) p n+2 - 2(n + 2)p n+1 + (n + l) 2 p n = — 

n + 2 

with the initial values po — 0, p± — 1. 

Proof. Applying Zeilberger's algorithm of creative telescoping [29, Ch.6] to the function 
F(n, t) we get for each n — 0, 1, 2, . . . the identity 

(26) F(ra + 2, i) - 2(n + 2)F(n + l,t) + (n+ l) 2 F(n, i) = G(n, t + 1) - G(ra, t), 
where 

t} = T(t+l)An r - n t + 3) ^ r(n ' = t{t2 " ( ^ + 3)t + + 2)) - 

To prove (26) it is sufficient to multiply both sides of (26) by r 2 (n-t+3)r(t+l)/r 2 (n+2) 
and after cancelation of gamma factors to verify the identity 

( n ^t + 2\ 2 

(n + 2) 2 -2(n + 2)(n-t + 2) 2 + (n-t + 2) 2 (n-t + l) 2 = - —!-r(n,t + l)-r(n,t). 

t ~\~ X 

Summing equality (26) over t — 0, 1, 2, . . . and taking into account that 

lim G(n, t) = 0, k — n + 3, n + 4, . . . , 

we get the difference equation for q n : 

q n+2 - 2(n + 2)q n+l + (n + l) 2 g n = -G(n, 0) = 0. 

In order to get the recurrence relation for the sequence p n , it is convenient to rewrite 
F n as an infinite sum 

t=0 
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taking into account that 

d 



lim ^-F(n, t) = for 
t^kdt 



k — n + 1, n + 2, 



Then differentiating (26) with respect to i and summing over t — 0, 1, 2, ... we get for 
each n = 0, 1, 2, 



(27) 
Since 



F n+2 - 2{n + 2)F n+1 + (n + l) 2 F n = lim G'(n, k + 1) - G'(n, 0). 



G'(n,f) = (n + l)! ; 



t 3 - (2n + 3)t 2 + tn(n + 2) 
r(t + l)r 2 (n-t + 3) 

3t 2 - 2t(2n + 3) + n(n + 2) 



(2^(n-* + 3) -^(f + 1)) 



we see that 



fe^oo 

and consequently, (27) becomes 



T(t + f)r 2 (n-t + 3) 
lim G'(n, k + 1) = and G"(n,0) 



n + 2' 



F n+2 -2(n + 2)F n+1 + (n+l) 2 F n 



n 



n + 2' 



n = 0,1,2,.... 



This implies that the sequence p n = F n + ^q n satisfies the same inhomogeneous recur- 
rence, and the lemma is proved. □ 

Now let us define a complex integral I n by means of the Meijer G-function: 



(28) 



t ^^0,2 ( n + l,n + l 
I n := n\ G 2a ( Q 



nr 
2tH 



T(t + 1) 



dt, 



where c > n is an arbitrary constant. 
Lemma 2. The following formula holds: 

F n = I n + o(\ 

where the constant in O is absolute. 
Proof. Since 

T 2 (t-n) 



as 



n — > oo, 



(29) 



r(* + i) 



V(t - n) 



T(t - n) 



by similar arguments as in the proof of Proposition 2, considering the integrand (29) 
on the rectangle contour with vertices c ± iN, —N — 1/2 ± iN, where N is a sufficiently 
large integer, we conclude that the integral (28) can be evaluated as a sum of residues 

11 



at the points n, n — 1, It is easily seen that the function (29) has double poles at 

the points 0, 1, 2, . . . , n and simple poles at —1,-2, Therefore, we have 



/ = „,« f res ( T ^- n A 

= »' 2 E 



fc=— oo 



t=k \sin nt F(n - t + 1) • t(l - t) ■ ■ ■ (n - t) 



+ Ul ( r (^+ l)r 2 (n-t + 1) (sinvrt) ) 



E 

fc=— oo 
n 



(n-k)\ ■ fc(l- fc) •••(«- fc) 



fc=0 



dt V.r(t + i)r 2 (n-t + i) 



I- 1^2 



'-l) k k\ 



- k (n + l) 2 ^(n + 2) 



Since 



E 

fc=0 



-l) fc fc! 



n + 2) 



^ fcj 
" ^ (n + 2) 2 



< 



x ^ 

^ fc! 

fc=0 



e, 



we get the desired assertion. 

Lemma 3. The following asymptotic formulas hold: 



F n = n\ 



q n = n 



n 1 / 4 \v £ 
e 2 ^ ( 1 



' n 1 ' 4 ^2^ 



- + 0(n- 1 / 2 ) 
0{n- l l 2 ) 



as n — > oo, 



as n — > oo. 



□ 



Proof. To prove the asymptotic formula for F n , we apply the saddle point method (see 
[14, Ch.6, Th.3.1]) to the complex integral (28). We begin by considering the quadratic 
polynomial 

(30) p(r) = n(r - l) 2 - r, 

which has two real roots 1 + ^ ± + 4^5, and denote the largest of them by 



to = 1 + 7T- 



2n V Ti 4n 



11 11 

1 + ^7TT + ^ 



n 



1/2 2 



n an 



3/2 



+ 0(n- 2 ). 



Then put c = nr and suppose that the contour of integration in (28) passes through 

the point utq. Then the asymptotic expansion of the gamma function (17) implies that 
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the following formula holds on the contour of integration: 

*r[t+i) = ^p{2io g r(t - n) - io g r(f + 1)} 

= exp |2 (* - n - ^) log(i - n )-(t+^ log(f + 1) 

+ 2n - t + 1 + log ^Phx + 0(rT 1/2 ) } . 

Making the change of variable t — nr and fixing the branches of logarithms that take 
real values on the interval (1, +oo) of the real axis, we get 

1 (t + 1) 

where 

/(t) = 2(r- l)log(r- 1) - rlogr-r + rlogn, #(r) = ^ _^ rl/2 - 
This gives 



(31) 4 = ^ / e^Wff(r) dr(l + O^)). 

^ J TO — ICO 

Now we determine the saddle points of the integrand (31), i.e., zeros of the derivative 
of the function f(r). It is easily seen that the zeros of the derivative 

/'(t) = 2 log(r — 1) — log r + log n 

are simultaneously the roots of the polynomial p(r) defined in (30). It is clear that 
t is the unique saddle point of the function /(r) and the contour of integration in 

(31) passes through it. Now for r = r + iy, — oo < y < +oo by the Cauchy-Riemann 
conditions, we have 

d d 

— Re/(r + iy) = -Im— /(r + iy) = -2 arg(r - 1) + arg r. 
Since Re To > 1, for y < we get 

and therefore, 

(32) Re/(r + iy) = (argr - arg(r - 1)) - arg(r - 1) > 0. 
dy 

This implies that Re/(r + iy) strictly increases as y increases from — oo to 0. If y > 
then we have 

7T 

< arg r < arg(r - 1) < - 

and hence 

— Re/(r + iy) = (argr - arg(r - 1)) - arg(r - 1) < 0. 

13 



7T 

- < arg(r- 1) < argr < 



Therefore, the function Re/(r + iy) strictly decreases as y increases from to +00. 
This proves that Re/(r + iy) attains its maximum on the vertical line Rer = r at the 
unique point r and we can apply the saddle-point method to calculate the asymptotics 
of I n , 

Jn = e^^^ir(r )|-to^)(l + 0^)). 

Since 

f"(r ) = ] + 1 u = 2n^(l + 0(n-^)), g(r ) = n 1/2 (l + 0{n^ 
t (t - 1) 

and 

/(r ) = -21og(r - 1) - r = logn - 1 - - i- + 0(n~ 3 / 2 ), 

we get 



_ Trv^n 1 / 4 



e nlogn-n-2VS( 1 + ( n -l/2)) = + Ofa" 1 / 2 )) . 

V e n i/4 



To compute the asymptotics of q n , we note that q n = n\C n (—l), where C n (x) = 
l.^ x n e -x^{n) _ Y^l =Q (%) is the Laguerre polynomial. Then the Perron asymp- 
totics for the Laguerre polynomials (see [40, p. 199]) yields 

and the lemma is proved. □ 

Theorem 1. Let {q n } n >o, {p n } n >o be defined by (23). Then q n G Z, D n p n G Z /or eac/i 
n — 0, 1, 2, . . . , and 

* - ™ = - ! *w G/i + ° ( "" /2) ) • = " ! 9 + 0(n " 1/2) ) 

as n 00. 

Corollary 1. TTie sequence p n /q n converges to Euler's constant sub- exponentially: 
^i- 7 = e" 4 ^(27r + 0(n- 1 / 2 )) as n ->• 00. 



4. A CONTINUED FRACTION FOR EULER'S CONSTANT 

It is interesting to mention on connection of our construction with a continued fraction 
of the Euler-Gompertz constant 

roc —x roo 

(33) 5:= / dx= / log(x + l)e~ x dx = 0.5963473623 .. . 

Jo x + 1 Jo 

found by Stieltjes in 1895 

1 00 / _ m 2 \ 1 12 2 2 3 2 

(34) 5 = - + K " ' 



2 m=i \2(m + 1) y 2-4-6 
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The n-th convergent of this continued fraction has the form where q n is defined in 
(23) and the sequence s n is generated by the homogeneous linear recurrence (24) with 
the initial conditions sq — 0, si = 1. Stieltjes derived (34) from the continued fraction 
of Gauss 

t-i / L I i . 1 \ - a(c-b) (H-l)(c-a+l) (q+l)(c-fe+l) 

/opr\ j (a,0+ 1;C+ ljZj = 1 c ( c +l) Z (c+l)(c+2) Z (c+2)(c+3) z 

1 } F(a,b;c;z) ~1- 1 - 1 - 1 

where F(a, 6; c; z) is the Gauss hypergeometric series 



rpt l \ (a)k(b)k k 

F ^ z ) = Y.^y^ z 



k=0 



We now recall how to deduce (34) following [42, Ch. 18]. Setting 6 = and replacing z 
by — cz in (35), and then letting c tend to infinity we get the expansion 

Here the equality sign can be considered as purely formal since the series on the left-hand 
side diverges excepting for z — 0. Although the power series is divergent, nevertheless 
the continued fraction converges to an analytic function of z, namely to 

i r° e~ u u a - 1 , 
q(z) = — . / era. 

yi ; r(a) y 1 + zu 

More exactly, one has [42, p. 353] 
1 f°° e" 

(36) — -*- = 7,-.7 



T\ajJ i + zu M ~T + T + T+ I +T+ l + 



which is valid for a G C, Re a > 0, and z G C \ (— oo, 0]. Replacing z by \jz and then 
dividing by z in (36) yields 

1 f°° e^u"- 1 , 1 a 1 a + 1 2 a + 2 

(37) — . / du = - - - - 

r(a) 7 + w z + l + 1 +z+ 1 + 

valid in the range of (36). Taking the even part of the continued fraction (corresponding 
to the 2n-th convergents) on the right-hand side of (37) gives 

(38) ' f£^4- 1 ° 2 < 1 + a » 3 ' 2 + a » 



T(a) J z + u z + a — z + a + 2 — z + a + 4 — z + a + G + 

valid in the range of (36). Now taking a = z = 1 in (38) leads to the expansion (34). If 
we take a = z = 1 in (36) we get the following continued fraction expansion: 

1 1 1 223344 

~T + T + T + T + T + T + T + T + T+"" 

Note that continued fraction expansion (36) can be found independently if we define 

1 f°° e~ u u a ~ 1 

hM: =W)L (TT^* 4 ' a>0 ' 
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and indicate the following relations, which can be easily deduced by integration by 
parts: 

h(a, b) = h(a, b + 1) + axh(a + 1, b + 1), 
h(a, b + l) = h(a + l,b + l) + bxh(a + 1,6 + 2). 
Then we can easily see that 

/(a,6+l)_l ax (b+l)x (a + l)x (b + 2)x (a + 2)x 

IT- T 



+ 



i + 



+ 



i + 



f(a,b) 

which for b = implies (36). 

On the other hand, we note that the convergence of the continued fraction (34) is 
slow enough to imply certain results on arithmetical nature of 5. As it was mentioned 
by Aptekarev in [6] 



n 

and therefore, one has only rational approximations to 5 : 

^-<5 = 0(e- 4 ^) as n ->• oo. 

q n 

Notice that the irrationality of the Euler-Gompertz constant 5 is still an open problem. 
All we know is that at least one of the two numbers 7, S must be irrational [6]. Indeed, 
from (33) we have 5 = e£'i(l), where 

t 



EAz) 



e 

— at, 
t 



arg^l < 7r, 



is the exponential integral. Using the formula (see [1, p. 228]) 
which implies that 



k=i 



k\k 



axgz\ < n, 



k=l 

(-l) fc 



(~l) fc 
k\k 



and the fact that the numbers e and Y^k=i ( ~k\l are algebraically independent (see [37, 
Ch. 7, Th. 1]), we get that the numbers e and S/e + 7 are algebraically independent 
and therefore at least one of the numbers 7, S must be irrational. Recently, Rivoal [36] 
gave a quantitative version of this statement by proving that for any e > there exists 
a positive constant d(e) such that for any integers p, q,r, q 7^ 0, one has 

> d(e) 
q ~ # 3 + £ ' 

where H = max(|p|, \q\, \r\). 

Using the following theorem from the general theory of continued fractions we will 
be able to find a continued fraction expansion (though not a simple one) for Euler's 
constant 7 whose n-th numerator and denominator coincide with the p n and q n , respec- 
tively. It seems to be the first non-trivial continued fraction expansion convergent to 
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V 

7 - - 

q 



r 

5- - 



Euler's constant sub-exponentially, the elements of which can be expressed as a general 
pattern. 

Theorem A. ([26, Th. 2.2]) Let {A n }, {B n } be sequences of complex numbers such that 

A_i = 1, A = b , B_i = 0, B = 1, 

and 

A nJ B n _! - 4,_iB n ^ 0, n = 0, 1, 2, . . . . 

T/ien t/iere exists a uniquely determined continued fraction bo + K(a„/6 n ) with nth 
numerator A n and denominator B n for all n. Moreover, 

b = A Q , ai = A 1 - A Q B U h = B u 

A n —\B n — A n B n _i A n B n -2 ~ A n ^2B n 

a n = — , b n = — , n — 2,3,4, 

A n -\B n -2 — A n -iB n -\ A n ^\B n -2 — A n -2B n -i 

Applying the above theorem to the sequences {p n }n>o and {q n } n >o we get the following. 
Theorem 2. Euler's constant 7 has the following continued-fraction expansion: 



1 1 5 36 15740 a. 



n 



7 £S a * n ' b * n ' 2 - 4 - 16 + 59 - 404 + " " " + b* n + " " " ' 

where 

a{ = l, a* = -l, a* = -5, a\ = 36, 0* = -15740, 
K 1 6* = 2, b* = 4, b* 3 = 16, b\ = 59, b* 5 = 404, 

and 

(n-l) 2 , A ,* 9A (ra-l)(ra-2) 
(40) < = - V - ; A n A n _ 2 , &; = n 2 A n _i + ^ ^ ^g n _ 2 , n > 6. 

i/ere g ra is a sequence of positive integers defined by (23) and A n is a sequence of integers 
generated by the third-order linear recurrence: 

, N (n - 1) (n - 2) A n+2 = (n - 2) (n + 1) (n 2 + 3n - 2) A n+1 
(41) 

-n 2 (2n 3 + n 2 - 7n - 4)A„ + (n - l)VA„_i, n > 3, 



wn£a £ne initial values Ai = —1, A 2 = —2, A 3 = —5, A 4 = 8. Moreover, A n is positive 
for any n > 4, and A 2ra 25 even /or any n > 1. 

Proof. Consider sequences {p n }n>o and {q n } n >o from (23) and put p_i = 1, g_i = 0. 
For n > 0, define n := p n ~iq n — p n q n -\. Then we have 

(42) = 1, Di=?)2 = -1, f 3 = -i 34 = 2, f 5 = ^- 

3 

From the recurrent equations (24), (25) we get the relation 

n — 2 

(43) t>„ = (n - l) 2 3 n _i H <? n _i, n > 1. 

77/ 

Since g n is positive for any n > and O4 > 0, it follows easily by induction that d n 
is positive for any n > 4. Taking into account (42), we get d n 7^ for n — 0, 1, 2, ... . 
Moreover, (43) implies that nfl n e Z for any n and nD n /2 G Z if n is even. 
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Now by Theorem A, we get that there exists a uniquely determined continued fraction 
b + K(a n /6 n ) with n-th numerator p n and denominator q n for all n, where b = p = 0, 
a>i=Pi- PoQi = I, h = q 2 = 1, and 

/ a a \ , Pn-2Qn ~ PnQn-2 . 

(44) a n = , b n = , n > 2, 

On-l n -l 

such that 

(45) 1 = K(a n /b n ) = - ^ ^ .... 
V ; 1 K n/ nJ 2 + b 2 + b 3 + 



From (43), (44) we have 



2 n - 2 g n _! 



a n = -(n-l) , n > 1. 

n n _i 

From the recurrent equations (24), (25) we obtain 

n — 2 

Pn-2<?n - PnQn-2 = 2nD„_ 1 H q n - 2 

n 

and hence 

(46) &n = 2n + — n>l. 

n D„_i 

Now let A n := nd n . Then we have A n e Z, A n is positive for any n > 4 and A 2n is 
even for any n > 1. Define po = Pi = P2 = 1, P3 = 3, P4 = 10, 

n(n - l)5„_i nA n _i 

Pn = ^ = 2 ' n 

and make the equivalence transformation of the fraction (45) by the rule (see [26, 
Th. 2.6]) 

a *n = PnPn-ldn, Ki = p n b n , 71=1,2,3,..., 

then using first several values of the sequence q n , 

g = l, q 1 = 2, q 2 = 7, g 3 = 34, q A = 209 

and formulas (44), (46) we get (39), (40). What is left is to show that the sequence A n 
satisfies the recurrence (41). From (43) we have 

A n - n(n - l)A n _! 

5n-i = o , n > 3. 

n — 2 

Substituting this expression in (24) we get the four-term recurrence relation (41), which 
completes the proof. □ 
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5. Concluding remarks 



In section 3 we considered the function F(n,t), which is a special case of a more 
general function (22), such that the sequences q n and F n defined by 

n oo n j oo , 

(47) g n = £ F(n, t) = J2 F(n, t), F n = £ ^F(n, t) = £ |f(u, f) 

t=0 t=0 t=0 t=0 

satisfy the second-order homogeneous and inhomogeneous linear recurrences (24) and 
(25), respectively. In this connection, it would be interesting to investigate the following 
question: does there exist a function F(n,t) of the form (22) such that both sequences 
(47) satisfy the same second-order homogeneous linear recurrence? If the answer is 
positive it would be therefore of interest to examine the convergence of — that could 
lead to a simpler continued fraction for Euler's constant. If this is not so, it would be 
curious to have a proof of this fact. Note that if we remove the condition 

s u 

£^£^- 

3=1 3=1 

in (22), the existence of such functions F(n,t) is known. For example, the rational 
function 

_ r(2t + 3n + l)r 6 (t + n)T 2 (t + 3n + 1) 
(n ' } ~ T(2t + 3n)r 6 (t + 2n+ i)r 2 (t) 

and the sequences 

(I 

q n = ^F(n,t), F n = 

t=i t=i 
that satisfy the second-order homogeneous linear recurrence 

(n + l) 5 F n+1 + 3(2n + 1) (3n 2 + 3n + 1) (15n 2 + 15n + 4)F n - 3n 3 (3n - 1) (3ra + l)F n _i = 
were used in [47] for constructing an Apery-like continued fraction for C(4). 
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